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^ ! Abstract 

The aim of this paper is to derive (by using two operators, representable by 
• a Jacobi matrix) a family of g-orthogonal polynomials, which turn to be dual 

, to alternative g-Charlier polynomials. A discrete orthogonality relation and a 

three-term recurrence relation for these dual polynomials are explicitly obtained. 
The completeness property of dual alternative q-Charlier polynomials is also es- 
tablished. 
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1. Introduction 
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It is well known that orthogonal polynomials are closely connected with spectral 
2 ! properties of symmetric operators, which can be represented in some basis by a Jacobi 



matrix, and with the classical moment problem. Namely, a spectrum of an operator, 
represented by a Jacobi matrix, is determined by an orthogonality measure for corre- 
^ ■ sponding orthogonal polynomials. If orthogonal polynomials admit many orthogonality 

relations, then the corresponding symmetric operator is not self-adjoint and it leaves 
room for infinitely many self-adjoint extensions. These extensions are determined by 
^ \ orthogonality measures for the appropriate orthogonal polynomials. 

I Contrary to orthogonal polynomials of the hypergeometric type (Wilson, Jacobi, 

Laguerre and so on), basic hypergeometric polynomials (or g-orthogonal polynomials) 
are not so deeply understood yet. These polynomials are collected in the g-analogue 
of the Askey-scheme of orthogonal polynomials (see, for example, [1]), which starts 
with Askey-Wilson polynomials and g-Racah polynomials, introduced in [2] and [3]. 
Importance of these polynomials is magnified by the fact that they are closely related 
to the theory of quantum groups. As an instance of such connection we refer to a paper 
[4], in which Al-Salam-Chihara g-orthogonal polynomials have been employed to con- 
struct locally compact quantum group SUq{l, 1). Another application of g-orthogonal 
polynomials is related to the theory of g-difference equations, which often surface in 
contemporary theoretical and mathematical physics. 

The purpose of present paper is to study completeness and duality properties of 
alternative g-Charlier polynomials. We shall show below that this originates a novel 
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type of g-orthogonal polynomials and a discrete orthogonality relation for them. To 
achieve this, we essentially use two operators Ii and q'^°, which are certain represen- 
tation operators for the quantum algebra [/^(sui^i) with a lowest weight (however, we 
do not use explicitly the theory of representations in what follows). The operator Ii is 
related to the three-term recurrence relation for alternative g-Charlier polynomials. We 
diagonalize the trace class operator Ii and obtain two bases in the Hilbert space: an 
initial basis and a basis of normalized eigenvectors of Ii . These bases are connected by 
an orthogonal matrix. The orthogonality relations for rows and columns of this matrix 
lead to orthogonality relations for alternative g-Charlier polynomials and for the func- 
tions, which are dual to these polynomials. We extract from the latter functions a dual 
set of polynomials and obtain a discrete orthogonality relation for them. As a result, 
one is led to the completeness property of dual alternative g'-Charlier polynomials. 
Observe that the present paper is a continuation of our research, initiated in [5] and 

[6]. 

Throughout the sequel we always assume that g is a fixed positive number such 
that q < 1. We use (without additional explanation) notations of the theory of special 
functions and the standard g-analysis (see, for example, [7] or [8]). 

2. Pair of operators (/i, J) 

Let be a separable complex Hilbert space with an orthonormal basis \n), n — 
0, 1, 2, • • •. We define on Ti two operators. The first one, denoted as q'^° , acts on the 
basis elements as q''°\n) — q^\n), and the second one, denoted as /i, is given by the 
formula 

Ii\n) = an\n+l) + a„-i|^ - 1) + K\n), (1) 



;i + ag2n+i)y'(i + ag2«)(l + ag2"+2) 
bn^q''\-——:j-7——^-^ + aq^ 



(1 + ag2n)(l + ag2n+i) (1 + ag2n-i)(i + ag2n) 

where a is any fixed positive number. Clearly, Ji is a symmetric operator. 

Since a„ — > and 6„ ^ when n ^ oo, the operator Ii is bounded. Therefore, 
we assume that it is defined on the whole Hilbert space Ii,. For this reason, Ii is a 
self-adjoint operator. Let us show that h is a trace class operator. For the coefficients 
an and hn from (1), we have an+i/dn — >■ q^^"^ and bn+i/bn q when n — * oo. Since 
< g < 1, for the sum of all matrix elements of the operator Ji in the basis |n), 
n — 0, 1,2, we have X^„(2a„ -|- 6„) < oo. This means that Ji is a trace class 
operator. Thus, a spectrum of Ii is discrete and has a single accumulation point at 0. 
Moreover, a spectrum of /i is simple, since Ii is representable by a Jacobi matrix with 
a„ ^ (sec [9], Chapter VII). 

To find eigenfunctions C,x of the operator Ji, IiC,x = X$,x, we set = Xln 
where /3„(A) are appropriate numerical coefficients. Acting by the operator Ji upon both 
sides of this relation, one derives that Yl'^=o Pn{X) {(in\n + 1) + On-il^T- — 1) + bn\n)) — 
Aj^^o /^n(A)|n), where a„ and 6„ are the same as in (1). Collecting in this identity 
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all factors, which multiply \n) with fixed n, one derives the recurrence relation for the 
coefficients /3n(A): 

/?„+i(A)a„ + /3„_i(A)a„_i + /3„(A)6„ = A/5„(A). 

The substitution 

reduces this relation to the following one 

-A„/3;+i(A) - C„/3;_i(A) + (A„ + C„)/3;(A) = A/3;(A), 

l + ^ _ 2n-l 



A — rJi ^ ' ""i n — nn^ 

"^n- q /-, , _ O.N/-, , _..9„.i\> — a? 



(l + a?2n)(i + a^2„+i)' " ^ (l + ag2n-i)(l + ag2n)- 
This is the recurrence relation for the alternative g-Charlier polynomials 

Kn{X;a;q) := 20i(g"", -ag"; 0; q,qX) 

(see, formulas (3.22.1) and (3.22.2) in [1]). Therefore, /3;(A) = K^{X;a;q) and 

For the eigenvectors we thus have the expression 

Since the spectrum of the operator Ii is discrete, only for a discrete set of values of A 
these vectors belong to the Hilbert space H. This discrete set of eigenvectors determines 
a spectrum of Ii. 

Now we look for a spectrum of the operator Ii and for a set of polynomials, dual to 
alternative g-Charlier polynomials. To this end we use the action of the operator 

J g--^o - a q'^° 

upon the eigenvectors which belong to the Hilbert space Ti. In order to find how 
this operator acts upon these vectors, one can use the g-difference equation 

(g-- - aq^)K^{X) = -aK^{qX) + X-^K^{X) - X-\l - X)K^{q-'X) (4) 

for the alternative g-Charher polynomials Kn(X) = Kn{X; a; q) (see formula (3.22.5) in 
[1]). Multiply both sides of (4) by d„ \n) and sum up over n, where dn are the coefficients 
of Kn{X; a; q) in the expression (2) for /3„(A). Taking into account the formula (3) and 
the fact that J\n) = (g~" — aq^)\n), one obtains the relation 

J6 = -a^,x + X-' 6 - A-^(l - A) e,-iA- (5) 
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We shall see in the next section that the spectrum of the operator /i consists of the 
points q^, n = 0, 1, 2, • • •. This means that J has the form of a Jacobi matrix in the 
basis of eigenvectors of I\] that is, the pair of the operators /i and J form a Leonard 
pair (see [9] for the corresponding definition). 

3. Spectrum of /i and orthogonality of alternative g-Chcirlier polynomials 

The aim of this section is to find, by using the operators /i and J, a basis in the 
Hilbert space 7i, which consists of eigenvectors of the operator 7i in a normalized form, 

and to derive explicitly the unitary matrix ?7, connecting this basis with the basis 
|n), n = 0, 1, 2, ■ ■ -, in Ti. This matrix leads directly to the orthogonality relation for 
alternative g-Charlier polynomials. For this purpose we first find a spectrum of /i. 

Let us analyze a form of the spectrum of /i from the point of view of the spectral 
theory of trace class operators. If A is a spectral point of the operator 7i, then (as it is 
easy to see from (5)) a successive action by the operator J upon the vector (eigenvector 
of Ji) .^A leads to the eigenvectors ^g^A, m = 0, ±1, ±2, ■ ■ ■. However, since Ji is a trace 
class operator, not all of these points may belong to the spectrum of Ji, since g^™A oo 
when m — > +oo if A 7^ 0. This means that the coefficient 1 — A' of iq-iy in (5) must 
vanish for some eigenvalue A'. Clearly, it vanishes when A' = 1. Moreover, this is the 
only possibility for the coefficient of ^g-iA' in (5) to vanish, that is, the point A = 1 is 
a spectral point for the operator Ji. Let us show that the corresponding eigenfunction 
^1 = belongs to the Hilbert space Ti. 

By formula (II.6) of Appendix II in [7], one has Kn{l; a; q) = 20i(9^", —aq^^; 0; q, q) = 
(— a)"g"' . Therefore, 

o_V" (-a;g)„(l +ag^") r^2(..^.._Y^ (-a; ?)„(! + ag^")a" 
^^1,^1; 2^ (l + a)(g;g)„a«g«(«+i)/2''"^''"'^^ ^ (1 + g)(g; g)„g-(3«-i)/2 " 

In order to calculate this sum, we take the limit d, e — > 00 in the equality 

(1 + aq^"){-a; q)n{d; q)n{e; q)n /aqy „(„_i)/2 ^ {-aq; q)oo{-aq/de; g)oo 
^^{l + a){-aq/d;q)n{-aq/e;q)n{q;q)n ^deJ {-aq/d; q)oo{-aq/e; q)oo 

(see formula in Exercise 2.12, Chapter 2 of [7]). Since 

hm {d; qUe; qUaq/de)^ = 

a,e— »oo 

we obtain from here that the sum in (6) is equal to {—aq; q)oo, that is, {^1,^1) < 00 and 
^1 belongs to the Hilbert space Ti. Thus, the point A = 1 does belong to the spectrum 
of the operator Ii. 

Let us find other spectral points of the operator Ii (recall that a spectrum of Ii is 
discrete). Setting A = 1 in (5), we see that the operator J transforms ^^o into a linear 
combination of the vectors and ^^o. Moreover, belongs to the Hilbert space H, 
since the series 

^ (l + a)(g;g)„a" 
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is majorized by the corresponding series (6) for ^^o. Therefore, belongs to the Hilbert 
space Ti. and the point q is an eigenvalue of the operator Ii. Similarly, setting X — q in 
(5), we find that is an eigenvector of Ii and the point q^ belongs to the spectrum of 
/i. Repeating this procedure, we find that all $,qn, n = 0, 1, 2, ■ ■ -, are eigenvectors of Ji 
and the set g", n = 0, 1, 2, ■ ■ ■, belongs to the spectrum of Ii. So far, we do not know 
yet whether other spectral points exist or not. 

The vectors ^^n, n = 0, 1,2, • • •, are linearly independent elements of the Hilbert 
space H (since they correspond to different eigenvalues of the self-adjoint operator /i). 
Suppose that values g", n = 0, 1,2, ■ ■ constitute a whole spectrum of Ii. Then the 
set of vectors ^^n, n = 0, 1, 2, ■ ■ ■, is a basis in the Hilbert space Ti. Introducing the 
notation := ^^k, k = 0, 1, 2, ■ ■ -, we find from (5) that 

J Sfe = -a Ek+i + q~^ Sfe - q-^{l - q^) Sfe_i. 

As we see, the matrix of the operator J in the basis S^, /c = 0, 1, 2, • • •, is not symmetric, 
although in the initial basis n — 0, 1, 2, • • •, it was symmetric. The reason is that 
the matrix (flmn) with entries a„m '■= f3rn{q^), Tn,n = 0, 1,2, where (3m{q^) are 
the coefficients (2) in the expansion = /5m(g")|^), is not unitary. This fact is 
equivalent to the statement that the basis S„ = ^^n, n = 0, 1, 2, • • •, is not normalized. 
To normalize it, one has to multiply S„ by corresponding numbers c„ (which are not 
known at this moment). Let S„ = c„S„, n = 0, 1, 2, • • •, be a normalized basis. Then 
the matrix of the operator J is symmetric in this basis. Since J has in the basis {S„} 
the form 

then its symmctricity means that c'^^Cna = c~^c„+ig~'^~^(l — q"'~^^), that is, Cn/cn-i = 
\/ aq^ / (1 — g"). Therefore, 

c„ = c(aV("+^^/V(?;?)n)^/^ 

where c is a constant. 
The expansions 

^qn{x) = Enix) = y^cA(g")|m) = ^a„„|m) (7) 

m m 

connect two orthonormal bases in the Hilbert space Tii. This means that the matrix 
{(imn), m,n = 0,1,2, ■■ ■, with entries 

n.n, ^ aV("+i)/2 (-a;g)^(l + ag^-) ^^^^ ^ , 
amn = cA{<l)=c(^ (g;g)„ {1 + a){q; qUa^qmim^m ) ^rni^ ■.a;q) (8) 

is unitary, provided that the constant c is appropriately chosen. In order to calculate 
this constant, we use the relation Y^^=o l^mnP = 1 for n = 0. Then this sum is a 

I/O 

multiple of the sum in (6) and, consequently, c = {—aq; q)oo ■ 
The matrix {cLmn) is real and orthogonal, that is. 
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Substituting into the first sum over n in (9) the expressions for amn, we obtain the 
identity 



oo 



n„n(n+l)/2 



(-ag™; 9)000™ {q;q) 



E 



Km{q'';a;q) Km'{q'']a]q) 



(1 + ag^"*) 



m 




(10) 



n=0 



which must yield the orthogonahty relation for alternative q'-Charlier polynomials. An 
only gap, which remains to be clarified, is the following. We have assumed that the 
points g", n = 0, 1, 2, • • •, exhaust the whole spectrum of /i. Let us show that this is 

the case. 

Recall that the self-adjoint operator Ji is represented by a Jacobi matrix in the 
basis n — 0, 1, 2, • • •. According to the theory of operators of such type (see, for 
example, [9], Chapter VII), eigenvectors S,\ of Ii are expanded into series in the basis 
|n), n = 0, 1,2, ■ ■ -, with coefficients, which are polynomials in A. These polynomials 
arc orthogonal with respect to some positive measure dfi{X) (moreover, for self-adjoint 
operators this measure is unique). The set (a subset of M), on which these polynomials 
are orthogonal, coincides with the spectrum of the operator under consideration and 
the spectrum is simple. 

We have found that the spectrum of Ii contains the points g", n = 0, 1,2, ■ ■ -. If 
the operator Ii would have other spectral points x, then on the left-hand side of (10) 
there would be other summands iixkKm{xk]CL]q) Kmi{xk]a]q), corresponding to these 
additional points. Let us show that these additional summands do not appear. To 
this end we set m = m' = in the relation (10) with the additional summands. Since 
^0(2^; a; g) = 1, we have the equality 



According to the formula for the g-exponential function Eq[a) (see formula (II. 2) of 



Appendix II in [7]), we have Xl^o '^"^"^"'^^''^^(^5 ^)n^ = {-aq;q)oc- Hence, Y.kf^^k = 



and all /^a,^ disappear. This means that additional summands do not appear in (10) 
and it does represent the orthogonality relation for alternative g-Charlier polynomials. 

As we have shown, the orthogonality relation for the alternative g-Charlier polyno- 
mials is given by formula (10). Due to this orthogonality, we arrive at the following 
statement: 

Proposition. The spectrum of the operator Ii coincides with the set of points g", 
n = 0,l,2,---. The spectrum is simple and has one accumulation point at 0. 

4. Dual alternative g-Charlier polynomials 

Now we consider the second identity in (9), which gives the orthogonality relation 
for the matrix elements amn, considered as functions of m. Up to multiplicative factors 
these functions coincide with the functions 




Fn{x;a\q) = 24>i{x, -a/x; 0; g, g"+ ), 



(11) 
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considered on the set a;G{g ™|m = 0,l,2,---}. Consequently, 

'^n^n(n+l)/2 (1 + ag^"^) ^ 



and the second identity in (9) gives the orthogonahty relation for a|g): 



m=0 



The functions Fn{x] a,b\q) can be represented in another form. Indeed, taking in 
the relation (III. 8) of Appendix III in [7] the hmit c — > oo, one derives the relation 

20i(g-"^,-ar; 0; g, = (-a)>-%0o(g-^ -ag'", - ; q,-qVa). 
Therefore, we have 

; a\q) = {-arq-'^Mq-"' , - ; -?7«)- (13) 

The basic hypergeometric function ^(po in (13) is a polynomial of degree n in the variable 
//(m) := Q'""* — ag"*, which represents a g-quadratic lattice; we denote it by 

d„(/x(m); a; q) := sMq'"", -« - ; 5, -?"/«) • (14) 

Then formula (12) yields the orthogonahty relation 

y il + aq^^)a^ q-^(^-^-^y^dMm))dr,Mm)) = (15) 

;^J-ag-;g)oo(g;g)m ^^^^ an^n(n+l)/2 V ^ 

for the polynomials (14) when a > 0. As far as wc know this orthogonality relation 
is new. We call the polynomials dn{fJ^{m);a;q) dual alternative q-Charlier polynomials. 
Thus, we proved the following theorem. 

Theorem. The polynomials dn{l^{iTi); a; q) , given by formula (14), orthogonal 
on the set of points fi{m) := q~"^ — aq"^, m = 0,1,2, ■■ ■, and the orthogonality relation 
is given by formula (15). 

Remark. The duality of polynomials is the well-known notion (see [7] and [11]) and, 
in particular, in the case of polynomials, orthogonal with respect to a finite number of 
discrete points, it reflects the simple fact that a finite-dimensional matrix, orthogonal 
by rows, is also orthogonal by its columns (c/ (9)). There is also an analytical way of 
deriving a dual orthogonality for polynomials, whose weight functions are supported on 
an infinite number of discrete points (see, for example, [12]- [14]). But this derivation 
is given in terms of a dual set of functions (for instance, their explicit form (11) for 
the case of alternative g-Charlicr polynomials is given above) and one needs to make 
one step further in order to extract an appropriate family of dual polynomials from 
these functions. Observe that in our approach to the duality of g-polynomials it is not 
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assumed that an initial family is orthogonal, this orthogonality property is straightfor- 
wardly derived. Besides, one naturally extracts from a dual set of orthogonal functions 
an appropriate dual family of ^-polynomials. 

Let be the Hilbert space of functions on the set m — 0, 1, 2, • • • with the scalar 
product 

where the weight function is taken from (15). The polynomials (14) are in one-to-one 
correspondence with the columns of the orthogonal matrix (Smn) and the orthogonality 
relation (15) is equivalent to the orthogonality of these columns. Due to (9) the columns 
of the matrix (a^n) form an orthonormal basis in the Hilbert space of sequences a = 
{a„ \n — 0, 1, 2, • • •} with the scalar product (a, a') = OnO^- This scalar product 
is equivalent to the scalar product (16) for the polynomials dni/J-im); a; q). For this 
reason, the set of polynomials dn{l^{Tn); a; q), n = 0,1,2, ■ ■ ■, form an orthogonal basis 
in the Hilbert space t^. This means that the point measure in (15) is extremal for the 
dual alternative q-Charlier polynomials dn{ii{rn);a; q) (for the definition of an extremal 
orthogonality measure see, for example, in [11]). This means that the dual alternative 
g-Charlier polynomials (14) form a complete system in the L^-space with respect to 
the point measure in (15). Observe that the completeness of alternative g-Charlier 
polynomials in the L^-space with respect to the point measure in (10) is a consequence 
of the fact that the operator (1) is bounded. 

A recurrence relation for the polynomials d„(/x(m); a; g') is derived from (4). It has 
the form 

(g-™-ag™)rf„(^(m)) = -ad„+i(/x(m)) + g^"4(/i(m)) - g-"(l - g")rf„_i(/i(m)), (17) 

where dn{^{m)) = dn{^Ji{m); a; q) . A g-difference equation for dn{^{m); a; q) can be 
obtained from the three-term recurrence relation for alternative g-Charlier polynomials. 
Note that for the polynomials (i„(/i(m); a; q~^) with g < 1 we have the expression 

ci„(//(m);a;g-i)=302(g-^-ag^g-"; 0,0; q,q). (18) 

However, the recurrence relation for these polynomials (which can be obtained from 
the relation (17), does not satisfy the positivity condition >l„C„+i > 0, that is, they 
are not orthogonal polynomials for a > (as it is the case for alternative g-Charlier 
polynomials). This positivity condition holds only if we demand a < 0. In this case, 
the polynomials (18) are the continuous big g-Hermite polynomials Hn{x;a\q) (for an 
explicit form of these polynomials see, for example, [1], formula (3.18.1)), which are 
orthogonal on a certain continuous set of points. 
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